Josephson junction transmission lines as tunable artificial crystals 
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We investigate one-dimensional Josephson junction arrays with generalized unit cells, beyond a 
single junction or SQUID, as a circuit approach to engineer band gaps. Within a specific frequency 
range, of the order of the single junction plasma frequency, the dispersion relation becomes gapped 
and the impedance becomes purely imaginary. We derive the parameter dependence of this gap and 
suggest designs to lower it to microwave frequencies. The gap can be tuned in a wide frequency 
range by applying external fiux, and persists in the presence of small imperfections. These arrays, 
which can be thought of as tunable artificial crystals, may find use in applications ranging from 
filters to the protection of quantum bits. 



I. INTRODUCTION 

The design of Josephson junction circuits in an appro- 
priate electromagnetic environment^ is currently of great 
interest in the context of superconducting qubits^. Such 
qubits are tunable artificial atoms. In analogy, one can 
envision long, periodic arrays^ of Josephson junctions, 
which are engineered metamaterials, or fully designable 
and tunable artificial crystals. In this context, we study 
the physics of regular, one-dimensional Josephson junc- 
tion arrays (Id-JJAs) with nontrivial unit cells, beyond a 
single junction or superconducting quantum interference 
device (SQUID). 

Regular JJAs have gained renewed interest: Id-JJAs 
have been the subject of study for the development of the 
10 Volt Josephson voltage standard^, and were more re- 
cently investigated for a current standard^, the transfer 
of quantum information^, and the readout and amplifi- 
cation of qubit signals^i^ at the quantum noise leveliSiii. 
2d- JJAs have long been of interest in the context of quan- 
tum phase transitions^^, and were recently investigated 
as negative permeability metamaterialsi^*^. 

Here we focus on the exploration of dispersion rela- 
tions and band gaps, well known from condensed matter 
and crystal physics. Very recently, crystal physics with 
Id-JJAs was studied in the regime where the Josephson 
junctions are treated as charge qubitsiS.. We consider the 
opposite limit of junctions in the phase regime, where 
junctions have a comparatively large area and therefore 
the advantage that they can be fabricated with a low rel- 
ative spread of parameters. Furthermore, we suggest two 
designs for which the band gap can be tuned by only one 
common applied magnetic field, which simplifies possible 
experiments. Junctions in the phase regime are approx- 
imately described by their linear behavior if the current 
flowing in the junctions is much less than the critical 
current. The array can then be regarded as a complex 
transmission line with a non-trivial, gapped dispersion 
relation. A resonator made with such a transmission line 
could find use in circuit cavity QEDiS,il2 for read-out of 
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FIG. 1: (Color online) (a) A model of a regular Id-JJA, where 
each unit cell j consists of a Josephson junction (modeled by 
tunnel element and capacitance) and a capacitance to ground, 
(b) A generalized model where each unit cell j consists of two 
different junctions and capacitances to ground. 



qubit signals. Nonlinear corrections to the model pre- 
sented here can be used to realize parametric amplifica- 
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and quantum noise squeezingi 



11.19.20 



We restrict ourselves to the linear regime, and focus on 
the possibilities of more complex unit cells and the pa- 
rameter regimes achievable in realistic experiments. We 
begin our studies with a unit cell consisting of two junc- 
tions (or SQUIDs) with different parameters, as shown 
in Fig. lb. In this case, we find a gap in the dispersion 
relation as shown in Fig. 2, i.e., a frequency region in 
which no propagating modes appear and the real part of 
the impedance vanishes. This gap appears at frequencies 
of the order of the plasma frequencies of the individual 
junctions. We investigate the parameter dependence of 
the gap and possibilities to lower it to the experimentally 
accessible frequency region, appropriate for qubit design. 
Then we discuss the connection between the real part 
of the impedance and the density of states for arbitrary 
unit cells in the linear regime, and present a more general 
treatment which no longer requires the unit cells to be 
identical. We use this general approach to show that the 
gap persists in the presence of a small parameter spread 
(5% standard deviation), and present simulations for a 
transmission experiment with realistic boundary condi- 
tions. 
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II. JOSEPHSON JUNCTION ARRAYS WITH 
TWO-JUNCTION UNIT CELLS 

A model of Id-JJAs in the linear regime is shown in 
Fig. la, where each unit cell consists of a Josephson junc- 
tion (with Josephson energy Ej and capacitance Cj ) , and 
a capacitance to ground Cq. In the linear phase regime, 
one can introduce the plasma frequency LUp — 1/ ^/LjCj, 
with linear Josephson inductance Lj — <i>§/(47r^i?j) and 
superconducting flux quantum <i>o — h/2e. Here, L] and 
Cj describe an intersite inductive and capacitive interac- 
tion, respectively, where we defined a site as the island 
between two junctions. The capacitance Co describes an 
onsite Coulomb interaction between Cooper pairs. We 
consider situations were resistive terms can be neglected. 
Also note that in the limit Cj ^ the model reduces to 
the discrete, lumped element model of a transmission line 
for transverse electromagnetic waves. We start our stud- 
ies with the generalization that the array consists of a lat- 
tice of unit cells each consisting of a basis of two Joseph- 
son junctions, for which we introduce an additional index 
1 or 2 to the parameters above, see Fig. lb. 

Let us first discuss a special case included in the gen- 
eralized model: If we choose Cqi = Cj2 = and regard 
Lj2 as a geometric inductance instead of a Josephson 
inductance, we recover a model containing one junction 
and an additional inductance in series. This model was 
studied earlier, where the inductance Lq was included 
in order to model the electromagnetic inductance of the 
JJA transmission line^^. This inductance led to a gap in 
the real part of the impedance, taken between input port 
and ground^ ^, which corresponds to a gap in the disper- 
sion relation as shown in Fig. 2. However, for typical 
parameters the gap appeared at approximately 10^^-10^'' 
Hz. While the lower frequency could, in principle, be re- 
duced by lowering the plasma frequency of the junction, 
the upper frequency extended beyond the range of valid- 
ity of the simple Josephson junction model used. With 
the two-junction model presented here, however, the up- 
per band edge can be reduced in frequency by orders of 
magnitude. 
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FIG. 2: (Color online) (a) Real part of the impedance (for an 
infinite array) and (b) dispersion relation. Both show a gap 
in the same frequency range. The impedance is normalized 
by Zji = ^Lji/Cji and the frequency by tJpi = 1/VLnCji. 
The length of a unit cell is called a, and k is the wave number 
of a traveling wave solution. Here we used parameters Cj2 = 
Coi ~ 0, where the model reduces to that of one junction 
and an additional inductance. Further, we chose parameters 
C02/CJ1 — Lj2/Lji = 0.1. 
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FIG. 3: (Color online) Dispersion relation for the model 
with two junctions when the junction parameters are (a) 
symmetric, L,j2/iji = 1, and (b) asymmetric, with weak 
{L32/L31 — 1.1) and stronger (I/j2/iji = 5) asymmetry. In 
all cases we used Cji = Cj2 and Coi/Cji = C02/CJ2 = 0.5. 



We shall now motivate the presence of this gap and 
then analyze its parameter dependence. First take an 
array with periodic boundary conditions and simple unit 
cells as in Fig. la, where each unit cell has only one de- 
gree of freedom due to loop constraints. Linearization 
of the equations of motion, in small value of phase dif- 
ference across each junction, approximates the system as 
coupled harmonic oscillators. A traveling wave ansatz 
yields a single branch in the dispersion relation. This 
branch will have an upper cutoff frequency due to the 
discreteness of the model but no second branch and no 
gap. If we now consider new unit cells consisting of two 
original cells, each of the new cells has two independent 
degrees of freedoms. This results in a representation of 
the dispersion relation where the original branch is mir- 
rored at half the Brillouin zone, and thus appears as two 
branches as shown in Fig. 3a. Here, we used the length 
a for the new unit cell. If an asymmetry of parameters 
is introduced within each unit cell, a splitting in the dis- 
persion relation occurs as shown in Fig. 3b. 

A geometric inductance Lq ^ introduces a much 
stronger asymmetry, which explains the wide gap and 
experimentally inaccessibly high frequency of the upper 
band edge (which tends to infinity for Lq/Lj — 0) in 
the model of Ref. [U With the two-junction model, we 
have a wide range of accessible parameters so that we 
can engineer the band gap in an appropriate frequency 
range. 

As shown in the appendix, the dispersion relation is 
given as 



B ^ y/B^ - AAC 
2A 2A 



(1) 



where 



A 



(Coi -1- Co2)(Cji -t- Cj2) -I- C01C02 + CjiCj2/3fc 



B = ^£L±^ + (^ + ^)/3, 



c 



with I/L12 = + l/ij2 and f3k 



(2) 

2[1 - cos(fca)]. 
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FIG. 4: (Color online) The upper and lower band edge in 
dependence of the ratio Lj2/Lji. We used Cj2/Cji = 0.5 
and Coi ~ C02 ~ 0.2Cji. Thin vertical lines mark where 
or ^2 change sign. 



Since A > 0, one sees from the defining Eq. ([T]) that 
0J^{k) > u)-{k) for any wave vector k. 

The lower and upper edge of the gap are defined as 



max{cL'_(fc)} , 

k 

min{cj+(fc)} , 

k 



(3) 



respectively. In the following we consider positive induc- 
tances and capacitances, for which the maximum of the 
lower band edge always appears at wave vectors ka = tt, 
where f3k = 4, while the upper band edge can appear 
at fca = or ka = tt, depending on the parameters, cf. 
Figs. 2b and 3b. More explicitly, we find 



OgV 



LO^{k — Ti / a) , 
uj+{k = 0) = 
a;+(fc — Tr/a) 



for 6-6<0 
for a-6>0 



where we defined Cs = C,n + Cj2 + CoiCo2/(Coi + C02), 
and 



= [Cq2Cj2 + C'oi(Co2 + C'j2)] Lj2 

— (Coi + C'o2)CjiLji 

= [C02CJ1 + C'oi(C'o2 + C'ji)] Lji 

— (C'oi + Co2)Cj2Lj2 ■ 



We show in Fig. [4] how the gap can be moved in fre- 
quency if one can control the ratio of effective Josephson 
inductances Lj2/iji. The gap vanishes, = Wgu, if 
both Coi = C02 = Co and 



L,nC,n — Lj2Cj2 -f- Co(ij2 — iji)/2 



(5) 



Kinks in the plot appear at the points, where either ^1 or 
^2 change sign. According to Eq. (|4]) this corresponds to 
switching the position of the minima of uj^{k) between 
k — and k = n/a, which can be shown to be realized by 
a flat (constant in k) upper branch i^+(fc) at these points. 



Tunable inductances as in Fig. |4] can be achieved by 
replacing one or both junctions in each unit cell by a 
SQUID, which consists of two parallel junctions with 
Josephson energy i?jo, and which is pierced by a mag- 
netic field B. The Josephson energy, which is inversely 
proportional to the Josephson inductance, can then be 
controlled as Ej = 2i!^jo cos(27r|_B|As/$o)- Here, As is 
the effective area of the SQUID with respect to the B- 
field. A design with one junction and one SQUID per 
unit cell has the advantage that one of the two induc- 
tances can be tuned continuously without changing the 
other inductance. A design with two SQUIDs per unit 
cell has different advantages. Clearly, one can then tune 
both Josephson energies. If one chooses different areas 
Asi and As2 for the two SQUIDs, a change in magnetic 
field AB = ^q/Asi leaves the Josephson energy of the 
first SQUID invariant, while it changes that of the sec- 
ond SQUID. In this sense, one can tune both Josephson 
energies independently with only one common magnetic 
field. Thus, this design is preferable in experiments which 
test many combinations of inductances (Lji,Lj2), while 
the aforementioned design is better if one needs a con- 
tinuous change of one inductance, which might become 
important for applications. 



III. MORE GENERAL UNIT CELLS 

Here, we will clarify why the real part of the impedance 
and the dispersion relation are connected as presented in 
Fig. 2, and show that this connection is of more general 
nature. At the same time, we will provide formulas which 
can also be used in the case of more general types of unit 
cells (e.g. with three or more different Josephson junc- 
tions or SQUIDs) and for the investigation of parameter 
spreads. 

If we allow an arbitrary combination of capacitances, 
inductances and Josephson junctions, the Lagrangian in 
the linear regime can always be written as 



L : 



].k=i 



(6) 



Here, M is the total number of independent variables, 
(4) and $j are independent flux variables. In the case 
of the two-junction unit cell we had M — 2N, with 
N the number of unit cells, and as variables we used 
the integrated voltage at the capacitances to ground, 
— Vj{t')dt'. Further, C is the capacitance matrix 
and is the inverse inductance matrix, which can con- 
tain both the kinetic inductance due to Josephson junc- 
tions and the geometric inductance. This is a problem of 
coupled harmonic oscillators, which can be diagonalized 
by the transformation $ ^ # = U^C^/^$. Here we 
took into account that the matrices C and L can always 
be chosen symmetric, and defined U as the matrix which 
has columns consisting of the normalized, real eigenvec- 
tors of the matrix fl'^ = C'^^^L-^C-^/^ . The trans- 



4 



formed Lagrangian is given as L = ^ X^jLiI^j ~ "^i^i]' 
where we introduced the eigenvalues of il^. These 
frequencies Wj resemble the dispersion relation for a reg- 
ular array, and can still be calculated in the presence of 
imperfections as shown in section IV. 

The equations of motion in this eigenbasis are decou- 
pled and given as = —u!^<bj. The Hamiltonian corre- 
sponding to the transformed Lagrangian is given as H = 

1 J2jiQ^ + '^i^i)' where Qj = $j is the conjugate vari- 
able to <&. For later convenience, this Hamiltonian can be 
rewritten in the standard form H = f^ji'^l^j + |) 

when creation and annihilation operators are defined by 
the equations 

Qj = *\/W2(a]-%) • (7) 

Now we proceed to consider the connection between 
the real part of the impedance, the density of states and 
the dispersion relation. For the frequency range at which 
the dispersion relation shows a gap, there are no un- 
damped, propagating modes and hence no transmission 
(in the limit of an infinite array). Then signals applied to 
the edge of the array having a frequency within this gap 
will be reflected back, and the impedance becomes purely 
imaginary. In general, the real part of the impedance can 
be expressed via the fluctuation dissipation theorem^^ as 

ReZH = ^(l-e-^''-)5yH . (8) 

Here, the spectral density is given as 5*^(0;) = 
I^oo s^^'^dt, where the voltage operator is 

given in the Heisenberg representation, V{t) = 
^iHt/hy^-iHt/h ^ assume a Boltzmann distribu- 

tion for the averaging, (O) = Tr{e-l^"0)/Z with Z = 
Tr{e~^^). Evaluating the trace in the energy eigenbasis 
H\n) — En\n) yields 

rn.n 

(9) 

In the case of a Hamiltonian consisting of M decou- 
pled harmonic oscillators discussed above, the states |n) 
and \m) are product states, and the energies En and Em 
depend on the occupation numbers of these states. For 
example, the energy of the state \n) — H^li Wj) is 

M ^ 

En = Y.^Ano^2^ ■ (10) 

i=i 

We shall assume sufficient spread of capacitances and in- 
ductances in an experiment, such that we do not have 
to consider degeneracies of the energies Cj = hujj. If we 
further consider for the operator of voltage in Eq. ^ the 



voltage at one site, Vk = ^'fc, we can rewrite the voltage 
operator as Vk = '^kj \/'2fiu!j{a^j—aj), with coefficients 

It is then straight forward to calculate the matrix ele- 
ments |(m |F| n)| in Eq. ([9]) to get 

j 

Using Eq. ([5]), the real part of the impedance is given as 

Re Zk{u; > 0) = ^ 6{iJ - e,/h) . (12) 

j 

The delta-function on the right side of the equation picks 
out only terms in the sum over j for which Cj = hw. The 
real part of the impedance is thus given as the number 
of states per frequency interval [in the continuum limit 
(X 1 / {du / dk)], multiplied by a weight factor \akj\^ which 
depends on the distribution of capacitances and induc- 
tances in the circuit. 



IV. EXPERIMENTAL CONSIDERATIONS 
A. Range of validity and experimental parameters 

The two branches in the dispersion relation and the as- 
sociated gap in frequency, where no propagating modes 
exist in the JJA transmission line, could be a useful prop- 
erty for the design of quantum circuits in the microwave 
region. The essential ingredient for realizing this gap is 
an inequality of the parameters for each of two junctions 
in the basis of the periodic structure, such that Eq. ([5]) 
is not fulfilled. Here we examine realistic designs, sub- 
ject to the constraints of fabrication, which can achieve 
this asymmetry. The designs naturally fall into two dif- 
ferent parameter regimes depending on the transmission 
line geometry used, coplanar or striplinc. 

When the JJA is made in a coplanar waveguide (CPW) 
geometry, where the ground plane is on the sides of the 
JJA, the parameter regime Cji, Cj2 ^ Cqi, C02 is easily 
realized. In this regime, the gap is vanishing or small 
when the plasma frequencies of the two junctions are the 
same, and one should thus aim for parameters Wpi ^ w-p2 
in order to observe the gap. When fabricating JJAs, typ- 
ically all junctions are made in the same process step, re- 
sulting in a tunnel barrier which is nearly uniform across 
the entire chip or wafer. In this case, the Josephson ca- 
pacitance Cj and inductance Lj will be proportional and 
inversely proportional to the junction area, respectively, 
and the plasma frequency ujp — \/ ^/LjCj will therefore 
be independent of the junction area. Thus, simply chang- 
ing the junction area in the fabrication process will not 
achieve the desired goal. 

Subject to the constraint of uniform tunnel barriers, 
there are two ways to bring down the plasma frequency. 
The first method is to increase Lj of one of the junctions 
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by forming a SQUID loop of this junction and applying 
an external magnetic flux. This method is attractive be- 
cause changing the external flux corresponds to tuning 
the frequency range of the transmission gap. However, 
with this approach, dropping the plasma frequency also 
drops the critical current of the transmission line, and 
therefore non-linear corrections will become important 
at much lower power. The second possibility to drop the 
plasma frequency is to fabricate an on-chip capacitance in 
parallel with each junction. This method will not cause 
a degradation of critical current, however, it does require 
more layers of lithography than the simple single layer 
process used in the shadow deposition technique. Fabri- 
cation with the Nb trilayer technique however provides 
this parallel capacitance naturally^^. 

When designing an array in CPW geometry, one 
finds that the characteristic impedance of the array 
is not well matched to the termination impedance. 
When the array is terminated with a direct connec- 
tion to an electrical lead, the termination of the ar- 
ray impedance at microwave frequencies will be approx- 
imately ZqI 2tt — 60ri, set by the free space impedance 
2'o = 377 ri. The transmission line impedance Zj^ of 
the JJA, which is the pure real impedance of an infinite 
array, is in the zero frequency limit given as .Za(O) = 

{Lji + Lj2)/(Coi + Co2)- It can be much larger than 
Zq in the CPW geometry, where Cq is relatively small, 
especially if Lj is made large by suppressing the critical 
current. We desire that Za ^ Rq = h/ie^ = 6.45kri, in 
order to avoid quantum fluctuations of the phase which 
are not included in our model based on classical phase 
dynamics. When Za 2> Rq, one finds that large quan- 
tum fluctuations of the phase result in a Coulomb block- 
ade, and our assumption of classical phase dynamics has 
completely broken down^^. 

An alternative route to circuit design is based on the 
stripline geometry, where the array is fabricated on top 
of a ground plane with a thin, insulating (non-tunneling) 
barrier separating the array islands from the ground 
plane. For the stripline geometry, one easily realizes the 
regime Cji,Cj2 ^ C'oi,Co2- In this regime, it is not 
necessary to have different plasma frequencies of the two 
junctions, and we find that a considerable gap in trans- 
mission also occurs if Cqi 7^ Co2- In this case the gap 
appears well below the plasma frequency. The condition 
of unequal capacitances to ground is easily realized in 
the stripline geometry, where the arrays are made as pla- 
nar circuits, with overlapping films forming the tunnel 
barriers. 

We have formulated a design for a JJA on an heav- 
ily oxidized Al ground plane, to be fabricated with the 
shadow evaporation technique. In our Al tunnel junc- 
tion fabrication, we find that it is possible to achieve 
plasma frequencies as low as a;pi/27r ~ LUp2/2T: = 33 
GHz. A design with large area base electrodes (2.5 /im x 
20 [iva) and long, narrow Dolan bridges (0.1 /im x 2.5 /im) 
with small overlap (0.1 /xm) after shadow evaporation 
can achieve the following parameters: Cqi = 0.68 pF, 



C02 = 3.4 fF, Cji = 2.0 pF, Cj2 = 12 fF, with the array 
critical current being dominated by the smaller junction 
2, Ic2 = 170 nA. For this design, we find that the lower 
gap edge comes down in frequency to Lo^i^jlix — 8.5 GHz, 
in a frequency range accessible to present day qubit de- 
signs or broad band transmission measurements. For this 
design, the transmission line impedance of the array is 
■Z^a(O) = 5317, which is well matched to the impedance 
of the input and output ports of an array with high fre- 
quency leads connected at each end. Such a design, with 
a rather low critical current and therefore strongly non- 
linear inductance, is ideal for the distributed parametric 
amplifier—. In the low power regime, where linear be- 
havior is expected, we find that such an array makes a 
good superconducting low pass filter, with a very sharp 
drop in transmission at 8.5 GHz in a design with only 20 
unit cells in series. 



B. Influence of parameter spread and transmission 

So far, we assumed that junctions can be fabricated 
identically if desired. However, in reality, there will be a 
spread of junction parameters in the fabrication. We ex- 
pect that the gap in the spectrum, impedance, and trans- 
mission will persist provided that the disorder is weak 
enough and/or the array is short enough that localiza- 
tion effects can be ignored. We investigated this further 
as shown in Fig. [5] (obtained by classical circuit theory), 
where we simulated a random spread of Josephson induc- 
tances with normal distribution and standard deviation 
5% in Fig. [5)3, d,f. Parameter spread, or disorder, vio- 
lates the condition of periodicity, and wave vectors are 
not well-defined. However, it is still possible to investi- 
gate the density of states, which is shown in Fig. [n^,b. 
Here we counted the number of states in a discretized 
frequency interval and considered 500 unit cells, as com- 
pared to 30 unit cells for Fig. [SJ:-f, in order to count a 
reasonable number of states. Despite the spread in pa- 
rameters, a gap can still be clearly observed in the density 
of states. 

In an experiment, it is easier to measure transmission 
than the dispersion itself. Furthermore, accurate bound- 
ary conditions on a finite length array become important 
for a real experiment. In Fig. [5j;,d we show the trans- 
mission for a finite array with boundary conditions de- 
fined by input and output leads with transmission line 
impedances Z^^ — Zout = SOfi. While this setup is not a 
cavity, it behaves like a resonator and supports standing 
waves because the input and output impedances Zout 
are not matched to the array impedance Za- The stand- 
ing waves can be seen by looking at the voltage at each 
site as in Fig.[S^,f. Note that because Za > Zin = Zout, 
the voltage antinode occurs in the middle of the array for 
the fundamental mode, opposite to standing waves in res- 
onators formed by a large point-like impedance at each 
end of a transmission line, where Zin, Zout > Za- Each 
standing wave condition is associated with a peak in the 
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FIG. 5: (Color online) (a),(b) Density of states. (c),(d) 
Transmission. (e),(f) Voltage at the sites. In parts (b), (d), 
and (f), a 5% standard deviation in Josephson inductances is 
used. Further, Lja = 0.25Lji, Coi = C02 = 0.2Ci = O.2C2, 
Zin = ^out, and Zin/ZA{0) = 0.2 . 



transmission as calculated in Fig. 6c, d for the case of 
no disorder, and 5% parameter spread, respectively. For 
frequencies inside the gap region, the transmission drops 
drastically which can be understood by an exponential 
decay of voltage amplitude from the edge of the array. 
When disorder or parameter spread is introduced, an in- 
vestigation of the eigenmodes shows localized states in 
the gap region, near the gap edge (not shown). However, 
the broad gap in transmission remains essentially unaf- 
fected by a 5% spread in parameters. Thus, the design 
of such a gap appears to be a robust and useful feature 
for quantum circuit engineering. 



V. SUMMARY AND OUTLOOK 

We theoretically investigated regular Josephson junc- 
tion arrays with generalized unit cells, e.g., unit cells con- 
sisting of two junctions with different parameters. In the 
linear approximation of the phase regime, we found a 
gap in the dispersion relation and in the real part of the 
impedance. This gap is not present in the linear regime 
of models with only one junction per unit cell, unless ad- 
ditional circuit elements in each unit cell are used. We 
derived the parameter dependence of the gap, and found 
that for a design with two different Josephson junctions, 
the gap appears at frequencies of the same order of mag- 
nitude as the plasma frequencies of the two junctions. 
We suggested how to lower these frequencies in an exper- 
imental setup in order to shift the gap to an accessible 
frequency range. By replacing one of the two junctions 
per unit cell with a SQUID, the gap can be tuned in situ. 
Furthermore, we derived the connection between the real 
part of the impedance, the density of states, and the dis- 
persion relation for more general unit cells, and found 
expressions which can be used when the periodicity con- 
dition is lifted. Using these expressions we showed that 
the gap persists upon a realistic parameter spread of the 
junctions. We also suggested a transmission experiment 
which we modeled with realistic boundary conditions. 

Our results could on the one hand be used for compar- 
atively simple demonstration of tunable artificial crys- 
tals with Josephson junctions, which is interesting from 
the point of view of fundamental physics. On the other 
hand, such tunable artificial crystals, embedded in su- 
perconducting circuits, could become important for new 
types of applications, for example, for specific types of 
frequency filters. 

One could also envision to protect qubits against de- 
coherence with such JJAs, by placing a Cooper pair box 
qubit in the middle of an array. If the energy splitting 
of the qubit lies inside the region of the gap, where no 
traveling modes are available, we expect the relaxation 
of the qubit to be strongly suppressed. Decoherence of 
a qubit is composed not only of the relaxation but also 
of the pure dephasing, where the latter time scale is typ- 
ically the critical, shorter one. However, recent experi- 
ments reached extremely high decoherence times, which, 
at least for part of the frequency range, appeared to be 
limited by the relaxations^. In this case, suppression of 
the relaxation with a properly engineered gap would al- 
low refined studies on remaining sources of decoherence. 

The model we presented in this paper is a linear anal- 
ysis of the JJA transmission line. The interesting ef- 
fects we describe arise due to the plasma resonance of 
the Josephson junctions when they are arranged in a 
discrete periodic structure. Many analogies exist be- 
tween the microwave phenomena presented here, and the 
emerging field of "plasmonics" , as represented by opti- 
cal plasma resonances in metallic and semiconducting 
nanostrucures^S.. However, the linear approximation is 
valid only when the currents fiowing in the junctions are 
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much less than the critical current. When this condi- 
tion is violated, nonlinear effects will appear, which can 
be very strong in comparison with dissipative effects in 
JJAs. These nonlinear effects give rise to a host of in- 
teresting phenomena, such as parametric amplification^. 
Here, the ability to match the JJA transmission line 
impedance for the case of the stripline geometry, with 
the electromagnetic transmission line impedance, leads 
to the possibility of a broad band parametric amplifier^*^. 
Another interesting nonlinear effect is the trapping or lo- 
calization of energy in discrete breather modes^l, where 
a gap in the dispersion relation is used to prevent the 
radiation damping of Josephson oscillations in a junction 
in the middle of the arrajs^^. We hope that the analysis 
presented here will aid the development of future exper- 
iments in these directions. 
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APPENDIX A: DERIVATION OF THE 
DISPERSION RELATION. 

The Lagrangian for the model consisting of unit cells 
with two types of Josephson junctions (Fig. [T)d) is given 



as 

N-l „ „ 
N-l „ „ 
AT-l 

+ E [-^Jl C0s{(t>j-i^2 - (t>j,l) + Ej2 COs((/)j,i - (t)j,2)] , 
3=0 

(Al) 

where we considered N unit cells, with periodic bound- 
ary conditions eliminated the Josephson phases (/^jj , and 
defined fluxes and phases $^,1/2 — 'i'o0j,i/2/27r at the 
capacitors to ground. 

From this Lagrangian, one can find the equations of 
motion, and, after the linear approximation, <^ 1, 
make a traveling wave ansatz. 

Here, we introduced a length a for the total unit cell, 
which results in a factor a/2 for a single junction. The 
equations of motion can be rewritten as a matrix M mul- 
tiplying the vector (u, v)'^ such that M(u, v)'^ = 0. Non- 
trivial solutions exist only when the determinant of M 
is zero, which results in the dispersion relation stated in 
Eq. ([T]) . The range of validity of the linear approximation 
is discussed in Section IV. A. 
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